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ABSTRACT 

We report on a deep, multiwavelength study of the galaxy cluster MACS J193 1.8-2634 us- 
| ing Chandra X-ray, Subaru optical, and VIA 1.4 GHz radio data. This cluster (z = 0.352) 

harbors one of the most X-ray luminous cool cores yet discovered, with an equivalent mass 
cooling rate within the central kpc is -700 M© yr" 1 . Unique features observed in the 
central core of MACS J193 1.8-2634 hint to a wealth of past activity that has greatly disrupted 
the original cool core. The X-ray and optical data suggest oscillatory motion of the cool core 
along a roughly north-south direction. We also observe a spiral of relatively cool, dense, X-ray 
\ emitting gas connected to the cool core, as well as highly elongated intracluster light (ICL) 

■ surrounding the cD galaxy. For a cluster with such a high nominal cooling rate, this cluster is 

missing the central metallicity peak almost always seen in cool core clusters, which suggests 
bulk transport of cool gas out to large distances from the center. Extended radio emission is 
observed surrounding the central AGN, elongated in the east-west direction, spatially coin- 
cident with X-ray cavities. The power input required to inflate these 'bubbles' is estimated 
from both the X-ray and radio emission to reside between Pj et ~4 - 14 xlO 45 erg s" 1 , putting 
it among the most powerful jets ever observed. This combination of a powerful AGN out- 
burst and bulk motion of the cool core have resulted in two X-ray bright ridges to form to 
the north and south of the central AGN at a distance of approximately 25 kpc. The north- 
ed ' ern ridge has spectral characteristics typical of cool cores: it contains low temperature, high 

density, metal rich gas and is consistent with being a remnant of the cool core after it was dis- 
rupted by the AGN and bulk motions. It is also the site of Ha filaments and young stars. The 
X-ray spectroscopic cooling rate associated with this ridge is -165 M yr _1 , which agrees 
with the estimate of the star formation rate from broad-band optical imaging (-170 M Q yr -1 ). 
MACS J193 1.8-2634 appears to harbor one of most profoundly disrupted low entropy cores 
observed in a cluster, and offers new insights into the survivability of cool cores in the context 
of hierarchical structure formation. 

Key words: X-rays: galaxies: clusters - galaxies: clusters: individual: MACS J193 1.8-2634 
- galaxies: clusters: general - galaxies: active - galaxies: cooling flows 
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1 INTRODUCTION 

* Email: sehlert@ stanford.edu 

f Einstein Fellow The relatively cool and dense gas at the centers of many galaxy 

$ Chandra/Einstein Fellow clusters emits copious X-ray radiation by thermal bremsstrahlung 
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and line emission. In the absence of external sources of heat- 
ing, this high emission should lead to very rapid cooling (t coo i < 
1 Gyr) and very high rates of mass deposition onto the cen- 
tral cluster galaxy (up to ~1000 M yr _1 ), in turn causing 
very high star fo rmation rates and stron g line emission around 
0.5-1.0 keV 



Table 1. Summary of the two Chandra observations of MACS J1931.8- 
2634. Exposure times are the net exposure after all cleaning and processing 
as described in Section [2] 



ir to rmation rates ana strong line emission a round 
(e.g. [Fabian & Nulsen|[l977l: ICowie & Binnevlll977l : 



IPeterson & Fabianll2006l : iMcNamara & Nulsenll2007h. The lack of 
such obvious observational signatures (e.g. Peterson et al.l 120011 , 
l2Q03h provides compelling evidence that some central source of 
heating must be present. The most plausible source of heating 
is feedback from the central Active Galactic Nucleus (AGN). 
Large X-ray cavities filled with radio plasma are clearly seen 
in many systems, which provide an expected so urce of heat- 
ing and turbulence (e.g. iBriiggen & Raised l2QQ2h . In systems 
such as the Perseus, V i rgo, Centaurus, and Hydra A Clusters 
(e.g. IFabian et all 120031: iForman et all )2005l: iNulsen et all [2005 : 



Fabian et all 120061: iForman et all 120071 : iSanders & Fabianl 12007 . 
20081 : ISimionescu et al.ll200^ Million et al. submitted), heating of 



the ICM has been argued to involve sound waves and weak shocks. 
Many questions remain as to the precise processes by which cen- 
tral AGN activity suppresses cooling flows, and the extent to which 
these feedback processes require finely tuned parameters. More 
ambiguous still are the nature and energetics of feedback mech- 
anisms at higher redshifts, as only a few systems with such high 
classi c al cooling ra t es have been studied i n detail JSchindler et al.l 
19971: lAllenl l2000l : iKitavama et aD 120041 : iBohringer et all 120051 : 



Qgrean et al.ll2010h . 



The galaxy cluster MACS J193 1.8-2634 is an extreme exam- 
ple of a cluster with a rapidly cooling core, making it an ideal sys- 
tem to test the limits of feedback mechanisms in galaxy clusters. In 
a short 12 ks Chand ra observation of MAC S J193 1.8-2634 taken 
in October of 2002 dAllen et al.ll2004l 120081) . X-ray cavities were 
detected surrounding a bright central AGN. The physical size of 
these cavities (~ 25 kpc) is similar to those observed in the nearby 
Perseus Cluster (IFabian et al.ll2003L l2006h . Indeed, with its lumi- 
nous cool core, visible central point source, and very large apparent 
cooling rate, MACS J193 1.8-2634 is in many ways a higher red- 
shift analog to the Perseus Cluster. Deeper observations of MACS 
J193 1.8-2634 were taken with Chandra in August of 2008, increas- 
ing the total clean exposure to ~100 ks. The combined X-ray data 
are presented for the first time here, and are complimented with 
optical (Subaru) and radio (VIA) observations. Our goals are to ac- 
quire a better understanding of the thermodynamic structure around 
the central AGN and the extent to which feedback from the central 
AGN counteracts this extreme cooling flow. 

The structure of this paper is as follows: Section [2 discusses 
the reduction of the Chandra data. Section [3] discusses the imag- 
ing analysis of the X-ray data, while Section 01 presents results on 
the mass profile and cooling flow of MACS J 193 1.8-2634. Sec- 
tion [5] discusses the spectral analysis. Sections [6] and [7] present the 
optical and radio data, respectively. Calculations of the energetics 
associated with the central AGN are discussed in Section [8j and 
the results are discussed and summarized in Section [9] Through- 
out this paper, a ACDM cosmology is assumed with £l A = 0.73, 
Q m = 0.27, and H = 7010ns- 1 Mpc" 1 . At the redshift of MACS 
J193 1.8-2634 (z = 0.352), 1 arcsec corresponds to 4.926 kpc. 



2 Chandra DATA REDUCTION AND PROCESSING 

Two Chandra observations of MACS J193 1.8-2634 were per- 
formed using the Advanced CCD Imaging Spectrometer (ACIS) in 
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Figure 1. Combined X-ray surface brightness profile for MACS J 193 1.8- 
2634 in the energy range 0.7-2.0 keV. The red points denote the total sur- 
face brightness, while the black points denote the background subtracted 
surface brightness. 



October 2002 and August 2008. The standard level- 1 event lists 
produced by the Chandra pipeline processing were reprocessed 
using the CIAO (version 4.1.2) software package, including the 
appropriate gain maps and calibration products (CALDB version 
4.1 .2). Bad pixels were removed and standard grade selections were 
applied to the event lists. Both observations were taken in VFAINT 
mode, and the additional information available in this mode was 
used to improve the rejection of cosmic ray events. The data were 
cleaned to remove periods of anomalously high background using 
the standard energy ranges and binning methods recommended by 
the Chandra X-ray Center. The net exposure times after processing 
are summarized in Table Q] 



3 X-RAY IMAGING ANALYSIS 
3.1 Surface Brightness Profiles 

Flat-fielded images were first created in the energy range from 0.7- 
2.0 keV for each of the Chandra observations. This energy range 
was chosen to minimize the impact of astrophysical and instrumen- 
tal background components. All subsequent imaging analysis was 
performed in this energy range. Surface brightness profiles were 
produced from these flat-fielded images centered on the central 
AGN (a(2000) = 19 /7 31 m 49*.6, (5(2000) = -26 /l 34 m 33 s .6). Each 
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surface brightness profile was fit with a King plus constant model 
in the range of 50 to 400 arcsec from the central AGN. The King 
plus constant model is explicitly given by the form 



3.3.2 Two -Dimensional Surface Brightness Modeling 



I(r) 



+ C 



(1) 



The best-fit constant for each observation was subtracted from 
each surface brightness profile, and Fig. Q] shows the exposure-time 
weighted average of the two surface brightness profiles. Over the 
radial range of 10-250 arcsec, the profile can be approximately 
described by a King model with r = 6.54 + 0.23 arcsec and 
P= 1.35 + 0.01. 

The background level determined from the surface brightness 
analysis described above was subtracted from each flat-fielded im- 
age on a pixel-by-pixel basis. These two images were reprojected 
into a common aspect solution and combined. This combined im- 
age, which is used in all subsequent imaging analysis, is shown in 
Fig.0 



3.2 Background Subtracted, Flat-Fielded Image 

On scales larger than -100 kpc, MACS J193 1.8-2634 exhibits 
a general elliptical symmetry, with the major axis aligned along 
the approximately north- south direction. There are asymmetries at 
these scales, however, that can be more readily observed by over- 
laying the X-ray surface brightness contours, seen in Fig. |2(b)| The 
isophotes are clearly not concentric, and the centroid appears to 
shift north and south of the central AGN with an amplitude as large 
as ~ 20-30 kpc. Inside the central 100 kpc, seen in Fig. |2(c)[ the 
morphology becomes considerably more complex. The bright cen- 
tral point source is surrounded by two bright ridges to the north 
and south. The northern ridge is brighter and has a relatively sharp 
boundary to the north, while the southern ridge is trailed by diffuse 
emission extending further to the south. These features, as well as 
the varying isophote centroids on large scales are indicative of past 
oscillations of the core, which currently appears to be moving to 
the north. 



3.3 Substructure Analysis 

3. 3. 1 Bandpass Filtering 

In order to better resolve small scale structures around the center 
of the cluster, a high-frequency bandpass filter was applied to the 
image shown in Fig. |2(a)| The functional form of the filter is given 
as 



F(k) = 



(if 



(2) 



We have set the scale length k = 10 arcsec. The transformed im- 
age is shown in Fig. |3(a)[ and shows the central AGN as well as 
the ridges to the north and south more clearly. Other features more 
apparent in this image are depressions in the X-ray brightness im- 
mediately to the east and west of the central AGN. Nearby sys- 
tems show clear cavities in the X -ray brightness near the central 
AGN, similar in shape to these (e.g. lFabian et al.l2003l ; lBirzan et al.l 
120041) . Unlike the cavities in those systems, however, there is no ev- 
idence in Fig. |3(a)| as to where the outer boundary of these cavities 
might lie. 



The X-ray image of Fig. |2(a)| was also fit with a two-dimensional 
elliptical beta model of the form 



1 + 



r(x,y) 2 
rl 



(3) 



where A is the central normalization, r is the core radius, and a is 
the power index. The term r(x, y) is given as 



r(x,y) 



[(x - xo) cos + (y- jo) sin 6] 2 
(y ~ yo) cos 6 - (x- xq) sin 6 



1 



which gives the distance of any point in the image (x, y) from a fixed 
center (x ,yo) with an ellipticity e. The position angle, 6, was also a 
free parameter in the fit. Since the image has already been back- 
ground subtracted, no further considerations for the background 
were included in the fit. The best fit parameters for this model are 
r = 3.21 +0.04 arcsec , £ = 0.290 + 0.004, 6 = 4.11+0.46 degrees 
and a = 1.139 ±0.004. 

The fit was then subtracted from the image, and the residu- 
als adaptively smoothed. The smoothed residuals image is shown 
in Fig. |3(b)| With the elliptical model subtracted, a spiral pattern 
beginning just east of the central AGN e merges. Such a pattern is 
expec ted to arise in off-axis mergers (e.g. Ascasibar & Mar kevitchl 
I2006L see also Section©. 



4 IMAGE DEPROJECTION ANALYSIS 

An image deprojection analysis of MACS J1931. 8-2634 wa s un- 
dert aken following the mann er described in lAllen et~al 1 d20Q8h (see 
also lSchmidt & AllerJl2007h . In brief, the azimuthally averaged X- 
ray surface brightness profile (centered on the central AGN) and 
the deprojected temperature profile are combined to simultaneously 
determine the X-ray emitting gas mass and total mass profiles of 
the cluster. We assu me that the dark+luminous mass distribution 
follows the model of iNavarro etaD i 19951 . 1 19971) (hereafter NFW 
profile) 



p(r) = 



p c (z)5 c 



(r/r s )(l + r/r s ) 2 



(4) 



where p(r) is the mass density, p c (z) = 3H(z) 2 /&7tG is the crit- 
ical density for closure at redshift z, r s is the scale radius, c 
is the concentration parameter (with c = r 2 oo/r s ) and 6 C = 
200c 3 / 3 [ln(l + c) - c/(l + c)]. The length scale r A is defined such 
that the enclosed mean mass density is A times the critical density 
of the Universe at the redshift of the cluster. For a given scale ra- 
dius and concentration parameter, a model temperature profile can 
be calculated and compared with the observed deprojected tem- 
perature profile. The observed deprojected temperature profile was 
azimuthally averaged and binned into annular regions with roughly 
5,000 counts in each (see Section O. The mass profile parameters 
were stepped over a range of values to determine best-fit values 
and uncertainties using a^ 2 minimization technique. This analysis 
assumed hydrostatic equilibrium and spherical symmetry, both of 
which are obviously not true in the vicinity of the central AGN. 
To account for these assumptions, the data for the central 50 kpc 
have been statistically down-weighted by adding 30% systematic 
uncertainties to the measured temperatures. 
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Figure 2. Combined, background subtracted, flat-fielded image of MACS J193 1.8-2634 in the energy range of 0.7-2.0 keV. a) The central 5'.7 x 5'. 7 of MACS 
J 193 1.8-2634. b) Same as in (a), but zoomed in by a factor of about 2 and overlaid with logarithmic surface brightness contours in blue. The centroids of these 
contours shift north and south of the central AGN by distances with an amplitude of up to -20-30 kpc. c) Same as in (a) but zoomed in by a factor of about 8 
to focus on the central AGN, the bright ridges about 25 kpc to the north and south, and the diffuse emission extending further to the south. 
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Figure 3. Images of MACS J 193 1.8-2634 that emphasize substructure, a) Image of MACS J 193 1.8-2634 after applying a high-frequency bandpass filter. This 
image better shows the central AGN and bright ridges to the north and south. Depressions in the X-ray emission also arise to the east and west of the central 
AGN. b) Image of MACS J 193 1.8-2634 after subtracting the best-fit elliptical-/? continuum model and adaptively smoothing the residuals. This image better 
shows a large scale (-200 kpc) spiral feature wrapping around the center of the cluster. 



Our best fit NFW profile has a scale radius of r s = 

Mpc, a concentration parame ter of c = 6.25^-^, 

and an equivalent velocity dispersion (Alle n etal.ll2002h cr' = 
^f50r s cH(z) = 966^g km s _1 . The integrated mass profile is shown 
in Fig. |4(a)| The enclosed mass within r 2 5oo = 505+3 9 kpc is 
M 25 oo = 2.64^ x 1O 14 M . 

The bolometric luminosity, cooling time, and equivalent mass 
deposition rate profiles have also been de t ermin ed using the calcu- 
lations described in detail in IWhite et al.ldl997h . These are shown 
in Figs. |4(b)] |4(c)[ and |4(d) | respectively. These profiles show that 
there is rapid cooling within the central 50 kpc. Within this region 
and in the absence of of balancing heat sources, a cooling flow 
with an equivalent mass deposition rate of M ~700 M© yr _1 would 
be expected. The bolometric luminosity within this region is ~1 
xlO 45 erg s _1 , and the cooling time is less than 1 Gyr. 



5 SPECTRAL ANALYSIS 
5.1 Methods 

The Chandra observations of MACS J193 1.8-2634 are sufficiently 
deep for high signal-to-noise spectra to be extracted from relatively 
small regions. This enables us to carry out detailed, spatially re- 
solved measurements of the thermodynamic quantiti es of the ICM . 
All spectral analysis was carried out using XSPEC dArnaudll2004i 
version 12.5). The backgrounds for all spectral analysis were ex- 
tracted directly from the science observations, specifically from a 
region roughly the same distance from the center of the detector 
as the cluster, but on the diagonally opposite chip. The background 
regions are devoid of point sources and cluster emission. 



5.1.1 Regions of Interest 



The spectral structure of MACS J 193 1.8-2634 was measured in 
both one- dimensional radial profiles and two-dimensional maps. In 
addition, the spectra of several regions of interest identified from 
the imaging analysis (specifically the central AGN and bright cen- 
tral northern and southern ridges) were also investigated in detail. 
To account for contamination by the central AGN, the 2.5 arcsec 
radius region surrounding the central AGN was excluded from the 
spectral profiles and maps. Region specific response matrices and 
ancillary response files were created for all spectra. 

Spatially resolved spectral maps for the cluster were extracted 
in regi ons determined by the contour binning algorithm of lSandersI 
(120061) . which creates bins of equivalent signal-to-noise, follow- 
ing contours in surface brightness. Regions were selected to have a 
signal-to-noise of 30, resulting in approximately 1,000 counts per 
bin. 

Azimuthally averaged spectral profiles were measured in an- 
nular regions with nearly equal numbers of counts. Initially radial 
profiles were made using annular bins that each contain roughly 
3,000 counts. Measuring the metallicity requires data with a higher 
signal-to-noise, so we also carried out a similar analysis using an- 
nular bins with roughly 10,000 counts each. The chosen center 
for all annuli was the location of the central AGN. The spectral 
properties of each annular region were measured both in projection 
and deprojection. Deprojection was implemented using the PROJCT 
mixing model in XSPEC. 
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Figure 4. Best fit profile from the image deprojection analysis, a) The best fit integrated gravitational mass profile in black, with the lcr confidence interval 
shown as the red curves, b) The integrated bolometric luminosity profile, c) The cooling time profile, d) The equivalent mass deposition rate profile. The high 
luminosity within the central 50 kpc leads to very short cooling times in this regi on and very large n ominal mass deposition rates. For comparative purposes, 
the nominal mass deposition rate in the Perseus Cluster is roug hly~400 yr'( Fabian et a l. 2002h. 



5.1.2 Modeling the Emission 



All spectral regions were initially modeled as a single tem- 
perature optically thin plasma using the MEKAL code of 
Kaastra & Mewel Jl993h incorporating the Fe-L calculations of 



Liedahl et al 



Jl995n~ and the phot o electri c absorption models of 



Balucinska-Church & McCammonl dl992h . W e used the determi- 
nation s of solar element abundances given bv lAnders & Grevessd 
dl989h . The abundances of metals (Z) were assumed to vary with 
a common ratio with respect to the Solar values. The single- 



temperature plasma model has three free parameters: the temper- 
ature (kT), the metallicity (Z), and normalization (K). 

In each region, the spectral anal ysis assumed a fixed Galactic 
absorption column of 8.3 x 10 cm 2 JKalberla et al.l2005h . consis- 
tent with the value measured direc tly from the X-ray spectra . The 
modified Cash statistic in XSPEC dCashlll979l : lArnaudlhoolh was 
minimized to determine the best fit model parameters and uncer- 
tainties. All uncertainties given are 68% (AC = 1) confidence in- 
tervals, unless otherwise noted. In the analysis of the bright ridges 
to the north and south of the AGN, the spectroscopic model in- 
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eluded an additional MKCFLOW component, appropriate for a sce- 
nario where gas is assumed to cool at constant pressure from an 
upper temperature down to a lower temperature. The upper tem- 
perature and metallicity of the MKCFLOW component were tied to 
their corresponding MEKAL components, and the low temperature 
was fixed to 0. 1 keV. 



5.1.3 Thermodynamic Quantities 

Several thermodynamic quantities can be calculated directly from 
the best fit MEKAL model parameters. The electron density (n e ), 
pressure (P), and entropy (S) of the ICM are derived from the 
MEKAL temperature (kT) and normalization (K) as 

9 4ttx 10 14 (1 +z) 2 D\K 

— (5) 



1.2V 



P. 

S - 



kTn Q 
: kTn: 2/3 



(6) 
(V) 



where the cosmological value of the angular diameter distance D A 
is 1016 Mpc at the cluster redshift. The volume of the region, V, 
is given in units of cm 3 . If the region had not been deprojected (in 
the case of the thermodynamic maps), the volume of the region was 
estimated as 



V : 



2 . 

mm 



(8) 



of any point in the region to the center of the cluster, respectively, 
and £ 2 is the solid- angle ext ent of the region in the sky jHenrv et all 
l2004l : lMahdavi et al.]l2005h . If the spectrum had been fully depro- 
jected (true for the annular profiles), the volume was calculated as 



(9) 



Since the uncertainties on the temperature are considerably larger 
than those on the density (> 10% for the temperature as compared 
to ~2-5% for the densities), the fractional uncertainties on the pres- 
sure and entropy given are similar to the corresponding temperature 
measurement. 



5.2 Spectral Results 

5. 2. 1 Thermodynamic Mapping 

The maps of temperature, density, pressure, and entropy in MACS 
J 193 1.8-2634 are shown in Fig. [5] One of the most notable fea- 
tures is the spiral of low temperature gas wrapping to the east 
and north of the central AGN seen in the large scale tempera- 
ture map. This is spatially coincident with the surface brightness 
excess seen in Fig. |3(b)| Such spirals can be due to either ram 
pressure stripping of an infalling subcluster core, or merger in- 
duced oscillatory motion of the cluster core. Both scenarios are 
seen in many nearby systems and hydrodynamic simulations (e.g. 
I Churazov eUdl 120031: lAscasibar & Markevitchll2006l: iDupke et al I 
20071: iLaeana et alJl20ld : lOwers et al J 120091 : IZuHone et alJbOld : 
iMillion et all201Gh . 

Zooming in on the temperature and entropy structure sur- 
rounding the AGN, Figs. |6(a)| and |6(b)| show that there are also 
two small regions of higher temperature to the east and west of 
the central AGN spatially coincident with the depressions in X-ray 
emission noted in Fig. |3(a)| Although the temperature of these re- 
gions is up to 2 keV higher than their surroundings, the statistical 



Table 2. Spectroscopic measurements of the X-ray bright ridges to the north 
and south of the central AGN. Temperatures are given in keV, metallicities 
are in solar units, and the cooling flow rate is given in M Q yr _1 . Upper limits 
are given at the 90% confidence level. 



Measurement 



North 



South 



kT 4.78 ± 0.64 



Z 0.53 + 0.11 0.22 ±0.09 

M 165+^ < 83 



significance is modest. The lowest temperature and entropy gas is 
located 25-30 kpc to the north of the AGN and spatially coinci- 
dent with the X-ray bright northern ridge. There is a similar ridge 
to the south which also has lower temperature and entropy gas, but 
at lower significance. 



5.2.2 Spectroscopic Cooling in the Bright Northern and 
Southern Ridges 

Since cooling flows are only expected in the inner, densest regions 
of the ICM, we have searched for the presence of spectroscopic 
cooling in the bright ridges to the north and south of the central 
AGN. The regions defined as the northern and southern ridges are 
shown in Fig. These regions were constructed to surround the 
brightest emission in the image shown in Fig. |2(c)| 

The two ridges are approximately equal distances from the 
central AGN, but have different spectral properties. The northern 
ridge has an emission weighted temperature of kT = 4.78 + 0.64 
and a metallicity Z = 0.53 + 0.1 1Z Q . This metallicity is higher 
than the average of this cluster of Z = 0.36 + 0.03Z©, discussed 
in more detail in Section [5.2.41 Adding a cooling flow component 
to the spectral model of the northern ridge, we measure a spectro- 
scopic mass deposition rate of M - 165 + 56 M yr _1 . Comparing 
the luminosity of the composite MEKAL+MKCFLOW model with 
the MKCFLOW component alone, the spectroscopic cooling of gas 
down to kT ~0.1 keV contributes ~30% of the emission from the 
northern ridge region. The southern bright ridge has a similar over- 
all temperature (emission weighted kT = 5.87^^), but a lower 
metallicity (Z = 0.22 + 0.09) and a 90% confidence upper limit 
on the spectral cooling rate of M < 83 M© yr _1 . These results are 
summarized in Table [2] 

We note that the MEKAL+MKCFLOW model offered a signifi- 
cant improvement to the spectral fit of the northern ridge (AC ~10 
with one additional fit parameter, a significant improvement at the 
~99.9% confidence level), almost identical to the AC improvement 
obtained with a two temperature (MEKAL+MEKAL ) fit, which in- 
troduces two additional fit parameters to the default MEKAL model. 
We also attempted to further resolve the spectroscopic cooling flow 
by using 2 or more MKCFLOW components with different tempera- 
ture ranges, but none of these more complicated models provided a 
further statistically significant improvement to the fit. For the south- 
ern ridge, the fit was not improved by the addition of either a cool- 
ing flow or a second MEKAL component. 
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Figure 5. Thermodynamic maps for MACS J 193 1.8-2634. a) The surface brightness image of Fig. |2(a)| with the field of view of the thermodynamic maps 
overlaid in blue, b) Temperature, kT, in units of keV. c) Pressure, in units of keV cm -3 , d) Entropy, in units of keV cm 2 . The lcr fractional uncertainties in 
the mapped quantities are < 20 per cent. The white region at the center of these maps is the exclusion region for the central AGN. There are a total of 60 
independent regions shown in these maps, which have a field of view approximately 200 arcsec (1 Mpc) in diameter. 



5.2.3 Azimuthally Averaged Thermodynamic Profiles 



The azimuthally averaged temperature profile for MACS J 193 1.8- 
2634 is shown in Figured] It is clear in this profile that the temper- 
ature does not decrease monoton ically towards the center as is ex- 
pected in a cool core cluster (e.g. I Allen et alj|200ll ; Vikhlinin et alj 



l2005h . but instead appears to increase in the central-most regions. 
This central jump in temperature is most clear in deprojection, 
where the innermost temperature is several keV higher than the 
measured temperature in the adjacent region. Such a sharp discon- 
tinuity in the temperature is often associated with shock heating, 
but the presence of the cooler ridges separated from the central 
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Temperature h_sokpc_ Entropy h_50kpc_ 

3.5 4.2 43 5.6 6.3 7 7.7 8.4 9.1 9.8 10 13 16 21 30 43 67 104 167 271 441 721 

(a) (b) 

Figure 6. The thermodynamic structure around the central AGN. These maps are identical to those shown in Fig. \5\ but zoomed in on the central AGN by 
approximately a factor of 4. a) Temperature, in keV. b) Entropy, in units of keV cm 2 . It is clear in both of these maps that the lowest temperature and entropy 
gas is located approximately 25 kpc to the north of the central AGN. 



AGN makes presence of shock heating ambiguous. To account for 
these cool ridge substructures, thermodynamic profiles were taken 
with the ridge regions excluded. These thermodynamic profiles are 
shown in Fig. [9] These profiles more accurately describe the av- 
erage thermodynamic structure of the ICM. After excluding the 
ridges, the deprojected temperature in the central-most region still 
does not decrease towards the center. The sharp temperature jump 
seen in the center of the deprojected temperature profile is less pro- 
nounced, but still present, and there are weak indications of dis- 
continuities in the central bin of the density profile as well. All of 
the thermodynamic profiles of Figure [9] are consistent with shock 
heating occurring in the central 20 kpc, but the small number of re- 
gions and relatively low signal-to-noise do not allow for us to claim 
an unambiguous detection. Such a jump in the central temperature 
could also be due to the presence of cavities devoid of ICM gas. In 
that case, the measured central temperature would be due to ICM 
gas in front of and behind the central AGN. This should lead to 
a flatter deprojected density profile in the central regions, a feature 
that is not observed. The innermost annular region of the thermody- 
namic profiles is also considerably larger than any apparent cavities 
in the X-ray emission, so it is unlikely that this temperature jump 
is due to projection effects involving cavities. 

Both the temperature and density profiles show other discon- 
tinuities within the central 100 kpc. The density profile decreases 
discontinuously at several locations while the temperature profile 
increases sharply at r ~70 kpc. These profiles lead to clear dis- 
continuities in the entropy and an unusually flat pressure profile. 
Such discontinuities arise from the presence of either cold fronts or 
weak (At < 2) shocks. Since the lower temperature gas on the inner 
edges of these fronts has a higher density, this profile is more con- 
sistent with the presence of cold fronts than shock heating. Heating 



from weak shocks cannot be ruled out, however, since the state of 
the ICM at earlier times before any bulk motion or shock heating 
is unknown. Deeper data are required to distinguish between these 
possibilities. 

5.2.4 Metallicity Profiles 

Surprisingly, the metallicity profile of MACS J193 1.8-2634 shown 
in Fig. [TOl exhibits no deviations from a constant metallicity of 
Z = 0.36Z© out to distances as large as 400 kpc. All attempts 
to fit the data to a multi-temperature MEKAL model did not lead 
to any significant increases in the central metallicity. A central 
metallicity pe ak is almost always observed in the profiles for cool 
core clusters (lAllen & Fa bian 199 8]; iDe Grandi & Molendi 200 ll: 
[ Leccardi & Molendil 120081 : iLeccardi et aD l2010l : lEhlert & Ulmerl 
120091) . The only region in the cluster that has a significant metal- 
licity enhancement is the bright ridge of X-ray gas to the north of 
the central AGN, where the metallicity is closer to Z = 0.5Z o . The 
unusual metallicity profile of MACS J193 1.8-2634 argues that the 
cool core may have undergone substantial stripping, likely associ- 
ated with its bulk motion. A dramatic example of stripping of a cool 
core due to bulk m otion has been rece ntly reported in the nearby 
Ophiuchus Cluster (iMillion et al.ll2010h . 



6 OPTICAL STRUCTURE OF THE CLUSTER CORE 

We observed MACS J193 1 .8-2634 in broad-band BVRIz filters with 
SuprimeCam on Subaru (Table [3]). The data were reduced with a 
dedicated weak l ensing and photome try pipeline based on the GA- 
BODS pipeline o f lErben et al.l(l2005l) as part of a larger cluster sam- 
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Table 3. Summary of the optical data. All images were taken with Suprime- 
Cam at the Subaru telescope. The quoted seeing values are those of the 
coadded images. 



2 6 14 25 39 56 76 99 126 155 

Figure 7. The regions defined as the bright northern and southern ridges, 
drawn over the image of Fig. [2(c)l in blue. These encapsulate the regions of 
the brightest X-ray emission in the cluster apart from the central AGN. The 
spectral properties of these regions are discussed in detail in Section [5221 
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Figure 8. The azimuthally averaged temperature profile of MACS J 193 1.8- 
2634. The red triangles denote the projected temperature profile, while the 
black squares denote the deprojected profile. The central region has a much 
higher temperature than the adjacent regions, which suggests that shock 
heating might be present. The radius of the temperature minimum corre- 
sponds to the distance of the northern and southern ridges. 



pie (von der Linden et al., in preparation). Figure fTTT a) shows a BRz 
three-color image of the central 2'.6 x 2'.6 of the cluster. 

Similar to the X-ray structure, the optical morphology of the 
cluster core exhibits clear structure in the North-South direction: 
the Brightest Cluster Galaxy (BCG) and the intracluster light (ICL) 
are highly elongated in this direction. In Fig. [TTT b) we highlight 
the extent of the ICL. The SuprimeCam images are very deep, and 
we have taken significant care in flat-fielding the images (von der 
Linden et al., in preparation). This allows us to trace the ICL to 
~28.3 mag arcsec -2 . Correcting for cosmological surface bright- 
ness dimming, this would correspond to 27 mag arcsec -2 at z = 0. 
North of the BCG, two stars mask the ICL, but there is evidence 
that the ICL extends slightly northwest of the stellar halos. To the 
south, the ICL appears to encompass the second brightest galaxy, 
located > 200 kpc south of the BCG. At these extremes, the de- 
tected ICL extends to ~ 200 kpc north and south of the BCG, 
but only to ~ 70 kpc east and west. The ICL of MACS J1931.8- 
2634 is thus highly elongated, with an axis ratio of ~ 0.3. At radii 
> 50 kpc, the ICL traces the overall gravitatio nal potential and is 
tied to the evolution history of the clus ter (e.g. iKelson et al.ll2002l : 
iNapolitano et al .1120031 ; iGonzalez et al.ll2005l) . The stars in the ICL 
are collisionless, like the dark matter, and thus the shape of the ICL 
should presumably reflect the shape of the core of the cluster dark 
matter halo. The elongation of the d ark matter distribution indicates 
the direction of the last merger fe.g. lRoettiger et al.ll997h « and thus 
MACS J193 1.8-2634 likely experienced a merger within the N-S 
direction. The timescale of this event is not clear, since the elon- 
gation persists for several Gyr dRoettiger et al.|[l997l ; iMoore et"aD 
120041) . A confirmation of this elongation of the dark matter halo by 
gravitational lensing is currently not possible - the ground-based 
data do not reveal any strong lensing features. Furthermore, the 
sightline to MACS J193 1.8-2634 is close to the galactic center 
(/ = 12.5669°, b = -20.09°), and the field is crowded with stars, 
which prohibit weak lensing shape measurements for an apprecia- 
ble number of background galaxies. 

Of particular interest is also the structure of the BCG. 
Fig. H3a) shows the central 30 arcsec x 30 arcsec of the Suprime- 
Cam BRz image. We see 'pink' and 'blue' filaments extending 
to the northwest and southeast, respectively. The 'pink' nebulos- 
ity to the northwest of the BCG signals emission in the blue (B) 
and red (z) bands. At the redshift of the cluster, the emission lines 
Ha, [Nil], and [Sll] are redshifted into the z band. The response 
of the CCD is not uniform across the filter due to the decreasing 
quantum efficiency at longer wavelengths, and Ha falls onto the 
most sensitive wavelength interval of the filter throughput. In order 
to single out the line emission, we subtract the adjacent band (/) 
from the z-band, scaling the /-band image such as to remove the 
large-scale light distribution from the dominant stellar population 
of the BCG. We identify the resulting emission as predominantly 
Ha line emission, although this needs to be confirmed spectroscop- 
ically. We also subtract the /-band from the B and /?-bands, in or- 
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der to better visualize emission not stemming from the underlying, 
old stellar population visible in the /-band. The resulting image is 
shown in Fig. [T2f b). The dominant feature here is a bright fila- 
ment to the northwest, of both Ha emission and blue emission. At 
the cluster redshift, the [Oil] doublet is shifted redward of the B 
filter, and thus the most likely interpretation of blue light is con- 
tinuum emission from young stars. To the southeast, a filament of 
bright blue emission is visible, which lacks strong Ha emission. 
Faint loops of Ha and blue light are also visible to the east and 
west. Overall, the structure is reminiscent of NGC 1275, the cen- 
tral galaxy of the Perseus cluster, which exhibits a rich system of 
Ha fi l aments, as well as filaments of young stars dConselice et al.l 
l200ll : iFabian et al.l Eo08). From the data at hand, however, it ap- 
pears that in MACS J193 1.8-2634, Ha filaments are accompanied 
by blue filaments (Ha alone would be visible as deep red, not as 
pink, in our images). The reverse is not always true, however: we 
do not detect Ha emission coincident with the blue filament to 
the southeast. For comparison, this is significantly different from 
NGC 1275 , where most Ha em ission is not accompanied by star 
formation dCanning et al]|2O10h . 

In Fig. \V2\ c) we overplot contours from the X-ray emission 
from Fig. [2b). We see that the brightest knots of the northwestern 
filament (which are bright both in Ha and blue light) coincide with 
the northern ridge identified in the X-ray imaging and thermody- 
namic mapping (Sect. 15. 2. TT ). The southwestern filament also could 
be associated with the southern ridge. The X-ray cavities, on the 
other hand, are at almost 90° angles to the bright filaments. 

In Fig. \l2\ d) we overplot contours from the radio emission 
(Sect. 13 . The brighter radio emission seems to coincide with the 
brightest Ha emission, whereas the fainter radio emission is elon- 
gated in an East- West direction, following the X-ray cavities. 

Finally, in Fig.Ete) and (f) we show the HST WFPC2 snap- 
shot exposure of MACS J1931.8-2634 (1200s, F606W filter). This 
allows us to see details at the very core of the BCG, which are un- 
resolved in the deeper SuprimeCam images. Apart from the central 
AGN point source, we see three bright knots to the north, along 
spiral- shaped filaments. 

Giant Ha filaments and on-going star formation, as evident 
by the presence of young, blue stars, are a common occurrence in 
the BCGs of cool-core clusters (Allenlll995l : [Crawford et alJfl999L 
lHatch et al.ll2007l : iRaffertv et al.ll2008ah . In the most Ha luminous 
syste ms, the Ha luminosity scales well with the star fo rmation 
rate (lAllenlll995l : ICrawford et al"1ll999l : ICT Pea et al ]|200i). How- 
ever, in general, star formation is not the only source of photo- 
ionizing radiation in BCGs, as is evidenced by Ha emission present 
even in the absence of young stars, the spatial offset of Ha fil- 
ament s and young star clusters observed in so me near-by sys- 
tems (ICrawford et al.l 120051 : ICanning et al.1l2010h . an d the higher 
[Nil] /Ha line ratios when young stars are not present jHatch et al.l 
120071) . In the BCG of MACS J193 1.8-2634, young stars are clearly 
associated with the Ha emission. The Ha emission is furthermore 
coincident with the "northern ridge" defined in Fig. [7] Assuming 
ratios of [Nll]/Ha=0.7 and [ Sll]/Hof=0.3 (typical of Ha-luminous 
BCGs, ICrawford et al.Hl999h , the Ha luminosity in this region is 
L(Ha) ~ 9 + 2 x lO 42 /*^ erg s" 1 . This wo uld make it the most 
Har-luminous BCG known. ICrawford et al.l furthermore find that 
the typical color excess in Ha luminous BCGs is E(B - V) ~ 0.3, 
albeit wi th large scatter. A ssuming Milky- Way type extinction, the 
standard Ikennicuttl |l998) conversion between L(Ha) and star for- 
mation r ate (which fits Ha luminous BCGs well, a gain with large 
scatter, ICrawford et aD Il999l : la Pea et al thus suggests 

SFR ~ 170M©yr _1 . This is in remarkable agreement with the X- 



ray mass deposition rate of the northern ridge (~ 165 + 60M o yr _1 ). 
However, one needs to keep in mind that the SFR value is only a 
rough estimate, using only broad-band imaging, and assuming av- 
erage values of quantities with large observed scatter. 



7 RADIO OBSERVATIONS WITH VIA 

1.4 GHz radio observations were made with the Very Large Ar- 
ray (VLA) of the National Radio Astronomy Observatory on 2006 
April 14. The data were obtained in A configuration, and the time 
on source was 54 minutes. A central radio source with flux density 
~70 mJy is clearly associated with the central AGN and the core 
of the elliptical galaxy. Approximately 45" to the south is a Nar- 
row Angle Tail (NAT) radio galaxy with flux density ~ 135 mJy. 
As shown in Fig. |13(a)[ the tails of the NAT are swept back to the 
south, in the same north- south orientation as the major axis of the 
central elliptical. The radio morphology of the central AGN (seen 
in Fig. 1 13(b)) is amorphous without the clearly defined jets or lobes 
that are found in many radio galaxies. Such amorphous radio struc- 
tures have been seen associate d with cP galaxies in cooling core 
cluste rs such as PKS074 5-191 feaum & O 'Pedll 99 ll : [Taylor et all 
1 19941) . 3C317 in A20 52 dZhao et al.ll 199 A and PKS 1246-410 in 
the Centaurus cluster (ITavlor et al.l20 02). Most likely the radio jets 
have been disrupted on small scales by dense gas. 



8 THE CENTRAL AGN: POWER AND ACCRETION 

The energy of the outburst occurring at the central AGN manifests 
itself as both radiative emission and as a jet that inflates cavities 
filled with radio plasma. The energy being input into both of these 
channels can be estimated from the observations. 



8.1 Estimating the Radiative Power 

The central AGN is sufficiently bright in X-rays to directly measure 
its spectrum and luminosity. Since its spectrum is presumed to be 
non-thermal in origin, we have modeled it with a variety of power 
law models. 

Source counts are extracted from a 2 arc sec region centered on 
the AGN. The spectral background is extracted from regions that 
are chosen to be near the AGN and avoid the bright ridges to the 
north and south, which are likely poor representations of the true 
background surrounding the AGN. The source and background re- 
gions are shown in Fig. [T4T a). The net spectrum for the AGN using 
these regions is shown in Fig [T4T b). Three different assumptions 
about the absorption have been examined: 1) that the absorption 
is due only to Galac tic contributions, and is fixed at the value of 
iKalberla et aD d2005h : 2) that the absorption is only due to Galac- 
tic contributions, but the column density nn is a free parameter in 
the fit; and 3) that the absorption includes both a fixed Galactic 
component and an intrinsic component at the redshift of the cluster 
with the intrinsic column density w H a free parameter. The results 
of these fits including the model fluxes and luminosities between 
0.7 and 8.0 keV are given in Tabled It is clear from the fit statis- 
tics that a larger absorption column density is favored (AC > 12, a 
significant improvement at a confidence level well above 99.9%). 
We conclude that the central AGN has a luminosity in the energy 
range of 0.7-8.0 keV of ~8 xlO 43 erg s -1 . Systematic uncertainties 
arising from the choice of the background region are estimated at 
-30%. 
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Figure 9. Azimuthally averaged thermodynamic profiles for MACS J193 1.8-2634 with the bright northern and southern ridges excluded, a) Projected (red 
triangles) and deprojected (black squares) temperature profiles, b) Deprojected density profile, c) Deprojected pressure profile, d) Deprojected entropy profile. 
There are clear discontinuities in both the temperature and density profiles at distances of r ~ 70 kpc and further discontinuities in the density profile, the 
origin of which may be either bulk motion of cold fronts or weak shock heating. 



8.2 Estimating the Jet Power 

The radio and X-ray emission both suggest that X-ray cavities 
filled with radio plasma reside to the east and west of the cen- 
tral AGN. The energy and power required to inflate and fill these 
cavities (4PV) can be estimat ed using the methods d escribed by 



Allen et all <2006h. (see a lso bunn & Fabianl 120041 : bunn et all 
2005: IChurazov et al.ll2002h . 



Neither the X-ray nor radio image provide an well-defined, 
unambiguous choice for cavity regions. The X-ray image shows 
two regions to the east and west of the central AGN that look like 
cavities, although the outer boundaries are difficult to determine, 
particularly given the presence of the bright ridges to the north and 



south. In nearby systems, depressions in the X-ray emission as- 
sociated with the full extent of the radio plasma-fille d cavities are 
typic ally visible only in very deep X-ray observations teirzan et al.1 
2008). The radio emission, on the other hand, has a much more nat- 
urally defined outer boundary, but our 1.4 GHz radio data for this 
z = 0.35 cluster do not show clear lobe or jet structure to the ra- 
dio source. The boundaries of this radio emission suggest that it is 
confined by the surrounding X-ray gas, but we caution that if the ra- 
dio plasma were 'leaking out beyond' the boundary of the cavities 
then estimating cavity volumes based on the radi o emission would 
lead t o overestimating the energy in the cavities dFinoguenov et al.1 
2008). With these uncertainties in mind, identical calculations were 
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Figure 10. Azimuthally averaged metallicity profiles of MACS J 193 1.8- 
2634. The black cur ve is the mean profile d erived for the Low Entropy Core 
(LEC) sample from lLeccardi et all ( 1201 Ol) scaled to the estimated ngo for 
MACS J 193 1.8-2634. The blue circles are the projected metallicity profile, 
while the red triangles are the deprojected metallicity profile. 



performed with two different sets of cavities: one set based on the 
structure of the radio emission (hereafter the radio cavities) and an- 
other set based on the apparent cavities in the X-ray images (here- 
after the minimal X-ray cavities). The true cavity volumes, 4PV 
enthalpy, and jet power are expected to lie somewhere between the 
values calculated from these two sets of cavities, both shown in Fig. 

Both the radio and minimal X-ray cavities were modeled as 
tri-axial ellipsoids with volume V = (4/3)7rrir w r d . The measured 
lengths r x and r w are the lengths of the axes in the image plane along 
and perpendicular to the jet axis, respectively. The final length, r d , 
is the axis of the cavity along the line of sight. The initial length 
of the axis along the line of sight was estimated as the smaller of 
the two planar axis lengths, but allowed to vary independently from 
them. Calculating the sound speed for the X-ray emitting gas with 
mean molecular weight \i = 0.62 and adiabatic index y = 5/3, 
we estimated the time scale for bubble formation as t age = (n/c s ). 
From this, we calculated the power required to inflate the cavities as 
roughly Pj et = (4PV/t age ). A Monte Carlo analysis was performed 
that drew the temperature, density, and all three spatial axes from 
independent Gaussian distributions. The assumed uncertainties on 
ri, r w , and r d were 20%, 30%, and 30%, respectively, leading to a 
systematic uncertainty of ~50% in the volume. Both the temper- 
ature and density were calculated from the major axis r x using a 
power-law parametrization of the profiles shown in Fig. [9] bet ween 
8 and 80 kpc (X(r) = X r a ,X e kT,n t ). From these variables the 
pressure, enthalpy, and jet power were then calculated. The prior 
assumptions and subsequent calculations of the cavity energetics 
are listed in Table [5] Although each 4PV calculation includes sys- 
tematic uncertainties of the particular cavity volume, the dominant 
uncertainty is identifying a particular set of cavities, which leads to 
an uncertainty in the 4PV enthalpy of approximately a factor of 8 
and jet power uncertian within a factor of roughly 4. 



We find that the total 4PV enthalpy between the two bub- 
bles, after accounting for all uncertainties, is approximately 1 - 
8 x 10 60 erg. This corresponds to a power input into the ICM from 
the jet of approximately Pj et ~ 4 - 14 xlO 45 erg s _1 . The mechan- 
ical energy going into inflating these cavities is nearly two orders 
of magnitude larger than the radiative emission of the AGN, and 
larger than the bolometric luminosity with in the central 50 kpc of 
the cluster. Based on the scaling relation of ICavagnolo et al.l j2010h 
and the radio luminosity, the inferred jet power is Pj et = 7.7 x 10 44 
erg s _1 , consistent with the lower range of jet powers measured here 
after accounting for the scatter. The temperature profile shows ev- 
idence for heating that goes out to roughly the same distance as 
the radio emission, so it is possible that the bubbles are in fact as 
large as the radio emission even though the X-ray depressions are 
much smaller in scale. The jet power derived from the radio cavities 
is comparable to the power input measured in the 200 kpc cavities 
of MS0735.6+7421 dMcNamara et al.l2005h . which was previously 
the system with the most powerful jets measured. 



9 DISCUSSION 

As a larger, more luminous, hi gher redshift analog to ne arby sys- 
tems like the Perseus Cluster (iFabian et allEool 120061) . MACS 
J193 1.8-2634 provides an extreme example of a cluster with a 
rapidly cooling core and powerful AGN feedback. This powerful 
AGN outburst combined with merger-induced motion has led to 
a cool core undergoing destruction to an extent previously unob- 
served in galaxy clusters. 

There is clear evidence that MACS J193 1.8-2634 has under- 
gone a merger event that induced large oscillatory motions of the 
core. On scales of r ~200 kpc, a spiral of cooler, denser gas seen 
in both the X-ray image and temperature map is observed to wrap 
around the core. Such spiral structures arise naturally from mergers 
and subsequent s loshing and are observable for several Gyr after 
the merger event jAscasibar & Markevitchll2006l) . Clear deviations 
from elliptical symmetry are seen in the isophotes, whose centroids 
shift with distance along the major axis to the north and south. On 
smaller scales near the core (r ~50 kpc), the X-ray data show that 
the core has a history of motion in the north- south direction. There 
are two bright ridges to the north and south of the central AGN. The 
northern ridge has a sharp northern edge (possibly a cold front), 
while the southern ridge has a diffuse tail of emission trailing its 
southern edge. Both of these features are consistent with the dens- 
est gas currently undergoing motion to the north. 

Our optical data independently suggest such a merger: the ICL 
around the cD galaxy is highly elongated in the north- south direc- 
tion. 5-band emission originating from a young stellar population 
is observed both to the north and south of the cD galaxy, while 
Ha emission originating from ongoing star formation is observed 
predominantly at positions coincident with the northern ridge. The 
presence of a young stellar population to the south without any Ha 
emission suggests that the primary region of star formation is mov- 
ing northward. 

In the midst of the motions of the cluster core, a powerful 
AGN outburst has taken place. The central AGN is bright in X- 
ray emission, with a luminosity in the energy band of 0.7-8.0 keV 
of ~8 xlO 43 ergs -1 . This AGN is surrounded by extended, amor- 
phous 1.4 GHz radio emission. The major axis of this radio emis- 
sion is spatially coincident with depressions in the X-ray emission. 
The physical extent of the cavities based on the observations is un- 
clear, but estimating cavities based on the X-ray and radio emission 
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Figure 11. (a): BRz image of the central 2'.6 x 2.6 (780kpc x 780kpc) of MACS J1931. 8-2634. The sightline to MACS J1931. 8-2634 is close to the galactic 
center, and thus most objects in this image are foreground stars. The box indicates the field of Fig.[l2] (b): The deepest band, R, smoothed with a (X/6 Gaussian 
kernel, and scaled to bring out low-surface brightness features. The black contour follows a surface brightness level of 28.3 mag arcsec -2 . Bright stars are 
marked as white circles. For the two stars north of the BCG, red circles indicate where the PSF surface brightness reaches -28.3 mag arcsec -2 . Note how the 
BCG envelope and the intra-cluster light are highly elongated in the N-S direction. 



Table 4. Spectral models for non-thermal emission of the central AGN. The three models are described in Section [ol The absorption column densities are 
given in units of 10 22 cm -2 and all fluxes and luminosities are given in the energy range of 0.7-8.0 keVand take into account their respective absorption. The 
units for the flux are 10 -13 erg cm -2 s -1 , and the units for luminosity are 10 43 erg s -1 . Errors listed are lcr confidence levels. The final column provides the 
C-statistic and degrees of freedom for the best fit model. 



Model Galactic Nu Intrinsic Nu Photon Index Normalization Flux Luminosity C/v 

Fixed Galactic Absorption 0.083 0.0 1.24 ±0.07 (2.05 ± 0.15) x 10" 5 1-76+g-} 4 7.27^ 8 44 519.81/497 

Free Galactic Absorption 0.42 + 0.11 0.0 1.70 ±0.16 (3.86+^ 4 ) x 10" 5 8.14+J-^ 505.68/496 

Free Intrinsic Absorption 0.083 °- 71 !o22 1.70 ±0.16 (3.83+g-||) x 10 -5 1.95+Jj- 2 ! 8.05+g$> 506.42/496 



gives a robust range for the 4PV enthalpy of the cavities and their 
corresponding jet power. The 4PV enthalpy of these cavities is suf- 
ficient to counteract the radiative losses from the central 50 kpc 
for 30-250 Myr. The jet power (~4 - 14 xlO 45 ergs" 1 ) identifies 
MACS J193 1.8-2634 among the most powerful cavity sources yet 
observed. The power input into inflating these cavities is approxi- 
mately 100 times larger than the measured radiative power of the 
central AGN, and four to ten times larger than the bolometric lumi- 
nosity of the cool core. 

Un like other more typical c ool co re clusters ( Allen & Fabianl 
19981: iDe Grandi & Molendil l200ll : IVikhlinin et all 120051 : 
Pratt et all 120071 : IWerner et al.l l2008h . the azimuthally- averaged 
metallicity profile for MACS J193 1.8-2634 shows no significant 
deviations from a constant value. Assuming that there once was 
a central metallicity peak in MACS J193 1.8-2634, this suggests 



that large masses of metal-rich gas has been stripped from the 
center of the cluster and displaced to the surrounding regions. 
The extent of transport required to account for the flat metallicity 
profile is strong evidence that the original cool core has undergone 
destructive stripping as it traversed from one side of the cD galaxy 
to the other. The only region with an exceptional metallicity is 
the northern central ridge. The central metallicity peak in cool 
core clusters is usually expecte d to be robust, even in clusters with 
powe rful central AGN activity teohringer et al.ll2004l ; lRasera et al.1 
l2008h . 

Merger and central AGN activity have resulted in the for- 
mation of two X-ray bright ridges roughly equal distances north 
and south of the central AGN. The X-ray bright ridge to the north 
of the central AGN has characteristics usually associated with a 
cool core. It is the location of the lowest temperature, highest den- 
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(a) (c) ' (e) 




Figure 12. Optical structure of the BCG of MACS J193 1.8-2634. (a): SuprimeCam BRz image of the central 30arcsec x 30arcsec. (b): For this image, the 
contribution from the old stellar population of the BCG (as traced by the SuprimeCam /-band image) was subtracted from each of the B,R and z images 
before combining them to a color image. This enhances the blue and pink features visible to the southeast and northwest of the central AGN. "Pink" signals 
contributions from predominantly the blue (B) and the red (z) channel. At the redshift of the cluster, the Ha line falls into the z-band, and thus this emission likely 
stems from Ha nebulosity surrounding MACS J 193 1.8-2634. The blue emission, on the other hand, likely signals a young stellar population. Interestingly, 
the Ha emission and young stars coincide in the northwestern region, whereas in the southeast Ha emission is absent, or significantly weaker, (c): Contours 
of the X-ray surface brightness map overlaid on the image in (b). The brightest knots in the northwestern filament coincide with the peak of the cluster X-ray 
emission north of the AGN point source. Thermodynamic mapping shows that this is also the coolest, densest part of the ICM, i.e. the cool core. A second 
peak (which also corresponds to cold, dense gas) is also seen close to the southwestern filament. The X-ray cavities, on the other hand, are located at ~90deg 
angles to the bright filaments, (d): Overlay of the radio emission on (b). (e): The HST WFPC2 snapshot, showing the same field of view as a), (f): The central 
7.5 arcsec x 7.5 arcsec of the HST snapshot. Note the bright knots in a spiral-like structure emanating to the Northwest of the central, brightest knot. 



sity gas in the cluster. It is also the location of the most metal 
rich gas. The X-ray spectrum of this northern ridge has a spec- 
troscopic cooling rate of M ~165 M©yr _1 down to 0.1 keV, in 
good agreement with the observed star formation rate estimated 
from the Ha emission, ~170 M©yr _1 . This Ha emission is ex- 
pected to be almost entirely from star formation, as it is roughly 
3 orders of ma gnitude larger than w hat expected from Case B 
Recombination djohnstone et al.|[l987l) . The spectroscopic cooling 
contributes a large fraction of the emission from this region, esti- 
mated at 30% based on the spectral modeling. The northern ridge 
might therefore be in the early stages of 'catastrop hic cooling' (e.g. 
iFabian & Nuisenlll977l : IPeterson & Fabianll2006h . The cooler ICM 
gas to the south of the cD galaxy has a lower metallicity and a cur- 
rent spectroscopic cooling rate consistent with zero, but is also the 
location of a young stellar population. This appears consistent with 
sloshing-induced stripping of the cool core throughout its oscilla- 
tions along the north- south axis. The asymmetric thermodynamic 
structure and different stellar populations of these ridges clearly 
indicate that there is no longer a single core of low entropy gas 



surrounding the cD galaxy. The extent to which AGN feedback 
contributed to the present-day thermodynamic structure of MACS 
J193 1.8-2634, apart from core sloshing, is beyond the scope of the 
current observations. The majority of the stripping and disruption 
of metal rich, low entropy gas from the original cool core could 
have been caused by the bulk motions, but it seems likely that 
the AGN outburst may have contributed to the separation of the 
preexisting core into two X-r ay bright ridges roughl y equal dis- 
tances from the central AGN (iGuo & Mathews! [20 lQh . This clus- 
ter's star formation is exceptional, especially as it does not satisfy 
al l of the empirical cond itions for central star formation discussed 
in iRaffertv et all <2008bl) . Although the cooling time and entr opy 
are below the thresholds put forward bv lRaffertv etafl d2008bl) . all 
of the systems with clear evidence for high star formation rates also 
have AGN jet powers smaller than the cooling luminosity, which is 
not the case in MACS J193 1.8-2634. 

Although the extent of stripping and cool core disruption 
in MACS J193 1.8-2634 is substantial, similar phe nomena have 
also been observed in the nearby Ophiuchus Cluster dMillion et al.l 
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Figure 13. 1.4 GHz radio emission in MACS J193 1.8-2634 observed with the VIA. a) X-ray image from Fig. [2(b)] overlaid with radio contours in magenta. 
The NAT galaxy to south of the center of the cluster is the source of the brightest radio emission in this cluster, b) X-ray image from Fig. |2(c)| with the radio 
contours overlaid in magenta. This figure shows the amorphous structure of the central radio source, which is clearly centered on the X-ray bright AGN. The 
radio contours are logarithmically spaced between 9 x 10 -5 and 0.1 1/y/beam. The beam size for this observation is 1.25 arcsec x 2.78 arcsec, and the position 
angle of the beam ellipse is 5.34 degrees. 





Energy (keV) 



Figure 14. Determining the radiative power emitted from the central AGN. a) The region of spectral extraction for the AGN (in blue) and background 
(magenta), b) The spectrum of the AGN with the best-fit power-law model including a fixed Galactic absorption and a free intrinsic absorption. The residuals 
of the fit are shown in the panel below. 



120101) . Many of the morphological structures seen in MACS 
J193 1.8-2634 are also seen in the Ophiuchus Cluster. Both clus- 
ters have clearly shifting isophotes, inner cold fronts, and comet- 
like diffuse emission trailing to one side of the inner cool core. The 
Ophiuchus Cluster also shows strong evidence for stripped core gas 
in the form of a metal rich ridge to the north of the cool core. These 



processes also appear to be occurring in MACS J193 1.8-2634, but 
with the added complications and energy of central AGN feedback. 

The extent to which a cool core can be disrupted or even 
destroyed by AGN feedback and merger induced oscillations has 
important implications for cosmologi cal studies with clusters (e.g. 
iBurns et al . 2008; Mantz et al.ll2010al lbh. Further observations with 
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Table 5. Calculations of the enthalpy an d jet power i n MAC S J 193 1.8-2634 for the cavities shown in Fig.[l5] These calculations follow the procedure described 
in Section [8!2l and are also discussed in I Allen et al.l {2006). The first seven rows are the prior assumptions going into the calculation while the final six rows 
are derived quantities. 



Parameter East (Radio) West (Radio) East(X-Ray) West (X-Ray) 

n(kpc) 20.2 ±4.0 30.0 ± 6.0 13.4 + 2.7 12.9 ± 2.6 

r w (kpc) 25.9 ±9.1 24.5 ± 7.3 9.4 ± 2.8 10.3 ±3.1 

r d (kpc) 20.2 ±6.1 24.5 ±7.3 9.4 ± 2.8 10.3 ±3.1 

£7o(keV) 5.30 ±0.25 5.30 ± 0.25 5.30 ± 0.25 5.30 ± 0.25 

a kT 

rc e ,o ( cm" 3 ) 3.47 ± 0.28 3.47 ± 0.28 3.47 ± 0.28 3.47 ± 0.28 

Qf ne -1.17 ±0.02 -1.17 ±0.02 -1.17 ±0.02 -1.17 ±0.02 

V (10 69 cm" 3 ) 1 - 25 -of L64 -0 65 °- n ± °- 08 °- 14 ± °- 08 

P(keVcm" 3 ) 0.50+0;^ 0.33^ 0.85 ± 0.20 0.81 ±0.25 

4PV(10 60 erg) 4.62 ±2.16 0.64 ± 0.32 0.85 ± 0.40 

cjkms" 1 ) 1170 ±30 1170 ±30 1170 ±30 1170 ±30 

*a ge (10 6 yr) 17.0 ±3.4 25.2 ± 5.02 11.4 ±2.1 10.7 ±2.1 

PjetClO^ergs- 1 ) 8.50 ±4.61 5.36 ± 3.06 1.90 ± 1.22 2.42 ±1.31 
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Figure 15. Estimating the jet power using the X-ray cavities and radio emis- 
sion. The image of Fig. |3(a)| is shown with the radio contours overlaid in 
magenta, the radio cavities overlaid in blue, and the minimal X-ray cavities 
overlaid in black. The determination of the 4PV enthalpy and jet power are 
derived from a Monte Carlo analysis that generates plausible elliptical cav- 
ities based on these ellipses. The true cavity volumes are expected to reside 
between the volumes calculated from these two sets of cavities. 



X-ray, optical, and radio instruments could provide many new in- 
sights into the extent that this cool core has been disrupted. Ra- 
dio observations at higher resolutions and lower frequencies could 
allow for a better understanding of the amorphous central radio 
source, in particular discerning the origin of the emission by mea- 
suring the radio spectrum. Higher resolution, deeper observations 
may also be able to resolve the jet and lobe structure of the central 
AGN source, which is critical for further constraining the magni- 
tude and origin of the AGN outburst. Optical spectroscopy would 
enable more precise measurements of the star formation rate in 
the different regions surrounding the central AGN, and perhaps al- 
low for measurements of the black hole mass. Optical spectroscopy 
would also allow for measurements of the emission line velocities 
within the cool core remnant, providing more details as to the ex- 
tent of the disruption of the cool core. Deeper observations with 
Chandra would provide better measurements of the unusual metal- 
licity structure and the extent that gas has been stripped from the 
cool core remnant. The distribution of regions that undergo cooling 
and the extent of that cooling could also be measured in more de- 
tail. Finally, a deeper X-ray exposure could shed new light on the 
thermodynamic structure of the cluster, in particular provide more 
compelling evidence for cold fronts and/or shock heating within the 
central 100 kpc. Simulations designed to reconstruct the thermody- 
namic structure of MACS J193 1.8-2634 may elucidate the nature 
of extreme sloshing and feedback, and perhaps also the future evo- 
lution of such a profoundly disrupted core. 
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